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Abstract 

The cosmic crystallography method of Lehoucq et al. P produces 



^ ■ sharp peaks in the distribution of distances between the images of 

I cosmic sources. But the method cannot be applied to universes with 

compact spatial sections of negative curvature. We apply to the these 
a second order crystallographic effect, as evidenced by statistical pa- 
rameters. 
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1 Introduction 



The method of cosmic crystallography was developed by Lehoucq, Lachieze- 
Rey, and Luminet P, and consists of plotting the distances between cosmic 
images of clusters of galaxies. In Euclidean spaces, we take the square of 
the distance between all pairs of images on a catalogue versus the frequency 
of occurence of each of these distances. In universes with Euclidean mul- 
tiply connected spatial sections, we have sharp peaks in a plot of distance 
distributions. 

It is usual to consider the Friedmann-Lemaitre- Robertson- Walker (FLRW) 
cosmological models of constant curvature {S^, E^, H^) with simply con- 
nected spatial sections. However, models with these spacetime metrics also 
admit, compact, orientable, multiply connected spatial sections, which are 
represented by quotient manifolds M = M/T, where M is S^, or H^, 
and r is a discrete group of isometries (or rigid motions) acting freely and 
properly discontinously on M. The manifold M is described by a fundamen- 
tal polyhedron (FP) in M, with faces pairwise identified through the action 
of the elements of F. So M is the universal covering space of M and is the 
union of all cells g(FP), g E T. 

The repeated images of a cosmic source is the basis of the cosmic cristal- 
lography method. The images in a multiply connected universe are connected 
by the elements g of F. The distances between images carry information 
about these isometries. These distances are of two types : type I pairs are 
of the form {g{x), g{y)}, where 

distance[g{x), g{y)] = distance[x,y], (1) 
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for all points x,y G M and all elements g G T; type II pairs of the form 
{x,g{x)} if 

distance[x, g{x)] = distance[y, g{y)], (2) 

for at least some points x,y & M and some elements g of F. 

The cosmic cristallography method puts in evidence type II pairs. These 
distances are due to Clifford translations, which are elements (? € F, such that 
Eq. Q holds for any two points x,y & M. Type II pairs give sharp peaks in 
distance distributions in Euclidean [HIS] and spherical spaces but they 
do not appear in hyperbolic space. This is illustrated in Fig. ^ior an FLRW 
model with total energy density Qq = 0.3 and having as spatial sections the 
Weeks manifold - coded m003(— 3, 1) in 5j and in Table I below - which is the 
closed, orientable hyperbolic manifold with the smallest volume (normalized 
to minus one curvature) known. The Bernui-Teixeira (B-T) function [B| is an 
analytical expression for a uniform density distribution in an open hyperbolic 
model. 

In hyperbolic spaces, the identity (or trivial motion) is the only Clifford 
translation. In this case, the cosmic cristallography method by itself cannot 
help us to detect the global topology of the universe. Several works have tried 
to identify multiply connected, or closed, hyperbolic universes by applying 
variants of the cosmic cristallographical method [3 |H1 El QQ] , most of which 
now rely on type I, in the absence of type II, isometrics. It is these variants 
that we call cosmic crystallography of the second degree. 

One of these 0, proposed by us, consisted of subtracting, from the distri- 
bution of distances between images in closed hyperbolic universes, the similar 
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— B-T function 
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Figure 1: Distribution of pair distances of cosmic images for the Bernui-Teixeira 
function, and for a universe with Weeks manifold as spatial sections and = 0.3. 

distribution for the open model. It did not pretend to be useful for the de- 
termination of a specific topology, but it might reinforce other studies that 
look for nontrivial topologies. 

Uzan, Lehoucq, and Luminet [9; invented the collect correlated pairs 
method, that collect type I pairs and plot them so as to produce one peak in 
function of the density parameters, for matter and Q\ for dark energy. 

Gomero et al. jTUj obtained topological signatures, by taking averages of 
distance distributions for a large number of simulated catalogues and sub- 
tracting from them averages of simulations for trivial topology. 

Here we introduce still another second order crystallographic method, in 
the absence of Clifford translations and sharp peaks. We look for signals of 
nontrivial topology in statistical parameters of their distance distributions. 

As commented above on Ref. [7], these methods are not as powerful as 
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the original Clifford crystallography, but will certainly be useful as added 
tools to help looking for the global shape of the universe. 



2 Simulations 



Let the metric of the Friedmann's open model be written as 



ds^ = a^{r]){d7]^-dX^) , 



(3) 



where a{7]) is the expansion factor or curvature radius, and 



dX'^ = dx^ + sinh^ xidO"^ + sin' 



^9d<p^) 



(4) 



is the standard metric of hyperbolic space H^. We assume a null cosmolog- 
ical quantity, hence the expressions for a{ri) and other quantities are as in 
Friedmann's open model - see, for example. Landau and Lifshitz 

To simulate our catalogues we assume for the cosmological density pa- 
rameter the values Qq = 0.3 and Qq = 0.4, with Bubble's constant Hq = 
65 kms~^Mpc~^. The present value of the curvature radius is ci(?7o) = 5513 
Mpc for Qo = 0.3 and 5954 Mpc for Qq = 0.4. 

To generate pseudorandom source distributions in the FP, we first change 
the coordinates to get a uniform density in coordinate space: 



with u{x) = (sinhxcoshx — x)/2 and v{6) = cos~^6. Our sources are then 
generated with equal probabilities in (m, v, 0) space, and their large scale 
distributions are spatially homogeneous. 



dV = sinh^ x sin d dx dO d(j) = du dv d(f) , 
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Table I. The eight manifolds used as spatial sections. Their names and 
data are taken from the SnapPea software. The last two columns give 
a number of cells that cover the observable universe. 



Manifold 


Number of cells 


Name 


Volume 




^out 


Qq = 0.3 


Qo = 0.4 


m003(-3,l) 


0.94 


0.52 


0.75 


747 


379 


m003(-2,3) 


0.98 


0.54 


0.75 


729 


357 


m017(-3,2) 


1.89 


0.64 


0.85 


463 


247 


m221(+3,l) 


2.83 


0.74 


0.94 


403 


199 


m342(-4,l) 


3.75 


0.77 


1.16 


451 


237 


s890(+3,2) 


4.69 


0.87 


1.38 


653 


273 


v2051(+3,2) 


4.69 


0.87 


1.38 


653 


273 


v2293(+3,2) 


4.69 


0.87 


1.38 


653 


273 



We did the simulations for eight spatially compact, hyperbolic models. 
Their space sections are the manifolds listed in Table I, which gives their 
names, volumes, and the circumscribing and inscribing radii {Xout Xin) 
of their FP's. These and other data on the manifolds were obtained from 
the SnapPea program 0. In the last two columns, we list the number of 
cells (replicas of FP) needed to assure a complete cover of up to a radius 
X = 2.335 for Qq = 0.3 and x = 1-996 for Qq = 0.4. This corresponds to the 
redshift {Z = 1300) of the last scattering surface. Aproximately 850 images 
were created in each simulated catalogue. Manifolds with different volumes 
will have different numbers of sources in their FP's. 
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Figure 2: The mean values for the eight studied manifolds, with Qq = 0.3. In this 
and remaining figures, a dashed line marks the value of the plotted parameter for 
the open model, and a solid line marks the Bernui-Teixeira value. 
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Figure 3: The mean values for our simulations with Qq = 0.4. 
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Figure 4: The standard deviations a for the simulations with Qq = 0.3. 
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Figure 5: The standard deviations for Qq = 0.4. 



8 



-0.54 p 
-0.545 - 
-0.55 - 

;/} 
(D 

1-0.555 - 

(D 
M 
00 

-0.56 - 



-0.565 - 



• 


mO03(-3,l) 


■ 


m003(-2,3) 


o 


m017(-3.2) 


o 


m221(+3,l) 


V 


m342(-4,l) 


> 


sS90(+3,2) 


A 


v2051{+3,2) 


< 


v2293(+3,2) 




open space 




B-T function 



2 3 4 

Manifold volume 



Figure 6: The skewnesses for VLq = 0.3. 
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Figure 7: The skewnesses for Qq = 0.4. 
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Figure 8: The peakednesses for Qq = 0.3. 
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Figure 9: The peakednesses for Qq = 0.4. 
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Ten different pseudorandom distributions of sources for each manifold 
were simulated. From the result of each simulation we calculated the distance 
distribution and then the latter 's mean value mv = E[X], standard deviation 
a = skewness^ sn = fi^/cr^, and peakedness^ pn = — 3, where 

X is the distance variable, E is the expected value operator, and /i^, = 
E[{X — mvY] is the //-moment about the mean value - cf. ^21, for example. 

For each of our chosen eight compact manifolds, we took the averages 
of these statistical parameters for ten simulations (obtained by varying the 
computer pseudorandom seed), to get their tendency in comparison with 
those for the simply connected case. 

We proceeded in two ways to obtain a distribution in the open model, 
for comparison with the compact cases. In one of them, a simulated cata- 
logue (real ones are not yet deep enough for use in cosmic crystallography) 
was obtained, with the same cosmological parameters and a pseudorandom 
distribution of sources inside the observable universe (redshift < 1300). In 
the other, the analytical Bernui-Teixeira function for a uniform distribution 

in simply connected universes |6j was used. 

^Skewness characterizes the degree of asymmetry of a distribution around its mean 
value. A negative value means a leftward asymmetry (as in Fig. 1), a positive value a 
rightward one. 

^Or kurtosis, which measures the relative peakedness or flatness of a given distribution, 
relative to the normal distribution. 



11 



3 Results and conclusions 



We compare the results for eight compact manifolds with those for the sim- 
ulated simply connected case, and with the B-T function. Figures 2-9 show 
these results for mv, a, sn, and pn, plotted vs. the manifold volumes, with 
the density parameters Qq = 0.3 and Qq = 0.4. It is of course premature 
to think of a functional dependence on volume (except for the obvious fact 
that, if the volume is so large as to enclose the observable universe, then the 
distribution of sources is indistinguishable from that for an open universe - 
see Fagundes for example). But we note a tendency of mv, a, and pn 
in compact universes to have values below those of the corresponding B-T 
function, and for sn to lie above that function. This could give us a signal 
of nontrivial topology in a future, realistic situation. 

We expect that the presence of type I pairs of distances is more evident 
for manifolds with smaller normalized volumes, or for models with smaller 
Qq. In fact our results are better for mv and a with Qq = 0.3, where the 
observable universe is bigger than for Qo = 0.4. This implies more copies 
of the FP, and hence more topological effects. The set of these statistical 
parameters may eventually provide a complementary indication as for the 
multiply connectedness of a possibly negatively curved cosmos. 
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